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Abstract 

The differential age data of astrophysical objects that have evolved passivelly during the history of 
the universe (e.g. red galaxies) allows to test theoretical cosmological models through the predicted 
Hubble function expressed in terms of the redshift z, H{z). We use the observational data for H(z) 
to test unified scenarios for dark matter and dark energy. Specifically, we focus our analysis on the 
Generalized Chaplygin Gas (GCG) and the viscous fluid (VF) models. For the GCG model, it is 
shown that the unified scenario for dark energy and dark matter requires some priors. For the VF 
model we obtain estimations for the free parameters that may be compared with further analysis 
mainly at perturbative level. 

Pacs: 98.80.Cq, 98.80.-k, 98.80.Bp 

1 Introduction 

The present cosmological observational estimations indicate that 95% of the matter in the universe is 
composed of exotic components which do not emit any kind of electromagnetic radiation in contrast 
to the ordinary forms of matter known in the laboratory. For this reason, these exotic components 
have been denoted the dark components of the universe. Ordinarily, these dark components are divided 
into two different ones: dark matter, exhibiting zero effective pressure, and directly connected with the 
formation of structure in the universe; dark energy, exhibiting negative pressure, responsible for the 
current stage of accelerated expansion of the universe. There are many candidates to represent dark 
matter (axions, neutralinos, etc.) and dark energy (cosmological constant, quintessence, etc.), all of 
them with important advantages and disadvantages. Since both components are detected only indirectly, 
through their gravitational effects, it is not excluded that the gravitational phenomena that arc interpreted 
as manifestation of dark matter and dark energy could also be explained by a deviation of the gravitational 
theory, i.e., general relativity. However, such proposal faces many theoretical and conceptual problems 
(conflict between local and global properties, etc.). In fact, all attempts to explain dark matter and dark 
energy may achieve remarkable success but, to our knowledge, all of them must also live with important 
drawbacks. For recent reviews, see references [TJ [H |31 HI [S] . 

In the last decade an interesting proposal to describe dark matter and dark energy has arised, the 
so-called unified model, where dark matter and dark energy are different manifestations of a unique fluid, 
keeping at same time general relativity framework untouched. The prototype of these unified models 
for dark energy and dark matter is the Chaplygin gas [6 , which has been generalized in many different 
ways [HI HI EH E] ■ if -essence models [P2] , which are based on modified expression for the scalar field 
kinetic term, can also unify the dark sector, and the Chaplygin gas is in fact a kind of iT-essence model. 
More recently, some special types of self-interacting scalar fields have been used to represent the unified 
scenario H3l H4l. 
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Dissipative fluid using the bulk viscous pressure |15j has also been considered as a possibility of 
implementing a unified model for the dark sector [16] , leading to very interesting results [171 118) . In the 
homogeneous and isotropic universe background a one-component viscous fluid shares the same dynamics 
as a GCG [19]. 

The unified models for the dark sector have been tested in many different ways, and the conclusion is 
that they can survive the observational tests, but in some cases specific priors must be introduced in order 
to avoid tensions among different observational tests. For example, confronting SNIa and LSS (Large 
Scale Structures) data, for the Generalized Chaplygin Gas (GCG) model, it seems that the unified model 
must be imposed from the beginning and the parameter a must be restricted to positive values, in order 
to alleviate the tensions in the parameter estimations. In particular, a must be restricted to positive 
values, in order to have a non imaginary sound velocity [20 , a crucial restriction when observational tests 
based on perturbative analysis are included. 

In this paper we will consider the H(z) observational data to investigate unified cosmologies using 
the fluid description for the cosmic components. Our candidates for the unified dark sector will be i) the 
Generalized Chaplygin gas, with equation of state 



where A and a are constants and p c is its density; ii) the bulk viscous fluid that is described by the 
Eckart's formula 



where O = 3H is the fluid expansion and £ = £(ft v ) = £,oPv ^ s the coefficient of bulk viscosity, written in 
terms of the density p v with constant parameters £o and v. 

We will not consider, initially, any other observational test, and the reason to concentrate on only one 
test is to verify to what extend there is potential tensions in the initial proposal of the model (for example, 
no need of an extra dark matter component). It is important, in order to have future reliable parameter 
estimations using all observational tests, to let a given observational test "to speak freely" , without further 
constraints. Later, as an application of our analysis, we will make a joint analysis using SN la, for the 
GCG model, for the best model obtained from the H(z) analysis. We use only one important restriction: 
the spatial section will be flat. The reason for this is twofold: the predictions for a flat universe are quite 
robust, and this seems to be corroborated by the position of the first accoustic peak in the spectrum of the 
anisotropy of the cosmic microwave background radiation |21j ; imposing a flat spatial section, we restrict 
the space of free parameters, avoiding unnecessary degeneracy as consequence of many free parameters. 
Our result at the end is very clear: the choice of the priors affects the final parameter estimation of any 
Chaplygin-based cosmology. In particular, concerning the unification program different priors over the 
parameter a can lead us to distinct conclusions. For the viscous model, we obtain estimations for the 
free parameters of the model, leading to results that are in qualitative agreement with the estimations 
obtained for the Chaplygin gas, but the unification scenario must be imposed from the beginning. 

The paper is organized as follows. In next section, we describe the models and the H(z) observational 
data. In section III, the GCG model and the viscous model are confronted against the H(z) data. In 
section IV, we present the final conclusions. 

2 The background dynamics of the unified models 

The GCG model is characterized by the equation of state ([1]) . The GCG is a realization of the energy- 
momentum tensor that appears in the right hand side of the Einsteins equations. The original interest 
for the Chaplygin gas came, from the theoretical side, to its connection to string theory [22 and, both 
from the theoretical and observational point of view, from the fact that it can unify dark matter and dark 
energy, representing both components at same time. Hence, it could seem unreasonable to consider the 
Chaplygin and an extra pressureless component. But, even in a unified scenario baryons must be present 
and, moreover, we would like to verify the priors and hypothesis that must be implemented in the model 
in order to have consistent predictions. 

If we add a matter component, which accounts for the baryons and an extra amount of pressureless 
dark matter, the flat Fricdmann's equation and the conservation law equations read, 
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The subscripts b, dm and c denote pressureless baryonic matter, pressureless dark matter and the Chap- 
lygin gas, respectively. In this sense, we begin with a three fluid model, where the total pressureless 
component is given by p m = Pb + Pdm- 

The conservation equations can be integrated leading to 
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where B is an integration constant. If we fix the scales such that ait = to) = 1, where to is the present 
time, pmo = pbo + PdmO 1S the present density for pressureless matter. Also, the GCG density can be 
re-written as 
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where p c o is the GCG density today, and A = A/ ' p)$ a is a parameter connected with the sound velocity 
of the fluid. 

The dimensionless density parameters are given by 
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where Hq is the Hubble parameter today. Hence, the Friedmann's equation may be written as, 
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Expressing now the scale factor in terms of the redshift z, a = 1/(1 + z), we find finally the expression 
for the function H(z): 



h(z) = H \ (n b o + n dm0 )(i + zf + {i-n 
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where i?o can be expressed in terms of the reduced Hubble parameter /i via Hq = lOQhKms Mpc 
and we have used that Q},o + Sl d mO + ^co = 1- The nucleosynthesis prior QbQ — 0.042 will be adopted. 

Equation ([9]) indicates how the Hubble parameter varies with the redshift. This is the quantity that 
we will be confronted with the observational data. 

In order to describe the background dynamics of the bulk viscous fluid we use the equation of state 
([2]) and we re-write equations ([3]) and (|4j) as 
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Remark that the continuity equation for the bulk viscous fluid depends on H . In the absence of pres- 

1/2 

sureless matter we have H ~ p v ' and the dynamics given by equations (|10H11|) corresponds to the set 
(j3]|4]) if we use the correspondence 

y=-{a+\). (12) 

Then, in a one fluid model, the GCG model and the bulk viscous fluid have the same background 
dynamics. 

Since pressureless matter was considered in the expansion factor (|10[) the correspondence described 
above does not occur within the cosmology considered here. Consequently, the background dynamics of 
the bulk viscous fluid will be different from the GCG model. The fractional density for the bulk viscous 
fluid will be given by the numerical solution of the continuity equation 
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where the new parameter £ can be expressed in terms of £0 via £ = — jj 2 ( j ■ 

Once we have solved the above differential equation for f2„ the final background dynamics is 

h v (z) = h [(n M + n dm0 ){i + zf + n v (z)] , (14) 

where the relation 51^0 + &dm0 + ^v0 = 1 still applies. 

3 Observational constraints using H(z). 

The Hubble parameter may be written in terms of the redshift parameter z as 

1 dz , . 

H= -. 15 

1 + z dt K ' 

Hence, to determine the function H(z) requires the knowledge of the derivative of z. This can be 
achieved by identifying some "clock" galaxies that exhibit a uniform distribution of star population. 
Special attention is given to red galaxies. This method has been proposed in [23] and employed to such 
galaxies in references [24, 25, 26]. We will use the data listed in reference [27] (see also [28]). 
The usual x 2 statistitics will be used, computing the quantity 

= (16) 

i ° l 

where the superscripts t and o denote the evaluated theoretical value and its corresponding observational 
result, while ai is the observational error bar. From this quantity, the probability distribution function 
PDF is constructed: 

V = Ae~x 2 / 2 , (17) 

where A is a normalization factor. 

The probability depends on the free parameters of the theoretical model. For the GCG model, these 
free parameters are h, Qdmo, ol and A. In our subsequent analysis, we will consider three main cases: 
the unification scenario, where the £X m o contribution is given just for the baryon component flbo] fixing 
Qdmo — 0.25 (the ACDM prior), where the GCG can be seen mainly as a dark energy candidate; leaving 
the parameter fldmo &s a free parameter. In this way we can consistently test the unification scenario of 
the GCG model. 

Choosing a prior is an important part of any Bayesian calculation. This includes to choose conveniently 
the interval for each free parameter. Concerning the parameters h, A and fldmo the natural choices are 

< /i < 1, < A < 1, < il dm0 < 0.958. (18) 

The background dynamics of the GCG model has a fundamental dependence on the parameter a. In 
principle, a can assume any arbitrary positive value. From the perturbative perspective, for example, 
a positive sound velocity requires a > 0. Even very large values, as a > 200, have been found in 
perturbative analysis 29, 20 . Thus, it seems that we do not need to fix an upper limit to a. On the 
other hand, the dynamics that has motivated GCG cosmologies, i.e., a matter dominated-like behaviour 
in the past followed by an accelerated one, occurs only if a > — 1. Also, some estimations in the range 
a < have been found using SNIa [30] ■ Hence, in order to explore different prior choices on a, we adopt 
the general priors a > a m i m , where a m i„ assumes four different values: -10.0, -2.0, -1.0, 0.0. 

We firstly discuss the results concerning the unification scenario (fldm — 0). Table 1 shows the one- 
dimensional estimations for each free parameter. The uncertainties are calculated at 2a. The results 
are not sensitive to the choice of the priors. There is a high concordance around the values h = 0.71 
and A = 0.97. Since the final estimation of any parameter depends only on the prior information of 
the remaining ones, in principle the estimated value for a does not depends on its own prior. This is a 
general result in bayesian analysis and is also valid for any other parameter. But this assertion is true 
only if the priors do not exclude regions of high probability for the parameter under evaluation. In any 
case, the uncertainties for the parameter a can be different as shown in table 1. This is due to the fact 
that the PDF is limited at a m i n - Consequently, the area below this function, which is used to calculate 
the uncertainties, depends on the prior. 
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Table 1: One-dimensional estimations of the parameters h, A and a for the unified GCG scenario (Sl^mo = 
0). 



Figured] shows a set of relevant PDFs obtained for the unified GCG model for the prior a m i n = —10.0. 
The one-dimmensional PDF for a (left panel) clearly shows a peak in its distribution at a = —0.20. 
The contours at 1,2 and 3 a of confidence level in the remmaining panels are obtained after the first 
marginalization, i.e., we integrate over A(h) in the central (right) panel. Since the contours a X h 
(center) and ax A (left) do not depend on the a-prior choice, we observe that a considerable probability 
is situated in the region a < 0. 
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Fi gure 1: PDFs for the unified GCG scenario (^d m o — 0) if a min — —10.0. 

In our next analysis we fix Vtdm — 0.25 (the ACDM prior). The final one dimensional estimations are 
displayed in table 2. The influence of the a-priors is weak and the estimations agree with h = 0.69 and 
A = 1. However, as shown in the PDF for a using a m i n — —10.0 (left panel of figure ^ there does not 
exist a peak in the a distribution. Its estimated value follows the a m i„ value what seems to contradict 
its ACDM limit i.e. a — 0. The two-dimensional PDFs in figure 2 show explicity that large negative 
values for a are preferred if the GCG is seen as a dark energy candidate. 
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Table 2: One-dimensional estimations of the parameters h, A and a for the GCG model with QdmO — 0.25. 



Leaving the parameter fldmo free to vary we have a consistent way to probe whether the unification 
scenario is favoured (fidmo — 0) or not (Qdmo > 0) . Table |3] shows that now the choice of the prior 
has a fundamental role. The unification scenario is favoured only if a m in = 0. Figure [3] shows the one 
dimensional PDFs for fi^mo f° r different values of a m i n ■ The two dimensional PDFs in figure _] show 
probability regions in the parameter space a x fldmo- When the marginalization process takes into account 
hypersurfaces corresponding to a < the probability for having higher values of fldmo increases. This is 
a strong evidence that the unification scenario is supported only if a m i n = 0. 
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Our second candidate for the unified dark sector is the bulk viscous model. For this model, we have 
still four free parameters when the spatial flat condition is imposed: Qdmo, h, £ and v. Assuming the 
unified condition i.e. Q v q = 1 — £lbo, which corresponds to the case (i) analysed above, the equation of state 
parameters are constrained as shown in figure [5J The estimations we have obtained are: h = 0.70^Q Qg, 
v = — 0.3ll^'g2 and £ = 2.43^o' 4 g, a t 2cr CL. We remark that the only difference between this model and 
the GCG one is due to the existence of pressureless matter. 

The unification scenario based on a bulk viscous model is well justified when we take into account 
the perturbative results [IS]. In fact, concerning the GCG model the matter power spectrum indicates 
that the unified scenario is not the preferred one (the estimated value for the pressureless matter is 
^dm ~ 1) t 20 . On the other hand, the same perturbative analysis for the bulk viscous models reveals the 
opposite result [18] . Hence, it seems contradictory to extend the analysis (ii) and (iii) to the bulk viscous 
model. Concerning the choice of the priors of the bulk viscous model there is no clear restriction from the 
perturbative dynamics. At first order the bulk viscous fluid behaves as a nonadiabatic component and its 
perturbations evolve in a very different way from the standard adiabatic scenario (which is the case of the 
GCG) |16j . Then, in principle, the parameter v can assume any real value. We remark, however, that at 
our knowledge no extensive quantitive parameter estimation has been made for the viscous unified model 
using matter power spectrum. 

From the analysis done so far using the H{z) data it seems that the unification scenario with the 
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Table 3: One-dimensional estimations of the parameters h, A and Qdmo- 




Figure 5: PDF for the bulk viscous case. 

prior a m i n = is the favoured one for Chap lygin- type cosmologies. Keeping this result in mind, we 
have also fixed the unification case for the bulk viscous fluid above. However, using other cosmological 
probes we can place stronger constraints on this unification case. In particular, a joint analysis using also 
the Sypernovae data will provide reliable constraints on the equation of state parameters a and A. In 
order to implement the joint analysis for the Generalized Chaplygin gas we rewrite x 2 (equation 17) as 
X 2 = Xsn + Xh' wnere Xh corresponds to equation (16). The Supernovae data set used in our analysis 
is the Union2 [31]. In figure|6]we assume the unification scenario (ftdmo = 0.0) and we show the contours 
for the equation of state parameters using only the Supernovae data (left panel) and the joint analysis 
(right panel). The contours obtained for the joint analysis arc wider as the errors added to the likelihood 
function (compared to the increase of the total % 2 ) become now larger 0. Also, note that the two samples 
predict a different estimated value for a yielding to a larger propagation of the uncertainties. For the 
joint analysis, we marginalize over the parameter A in order to obtain an estimation for a and we obtain 
a = O.lOig 57 which agrees (including the errors) with the result considering only the H(z) data. 

One important issue concerning parameter estimations for unified models is to discuss the constraints 
from the CMB physics. The allowed parameters we have found in this work is an acceptable fit CMB 
data according to previous estimations 32, 33 , which predict a ~ 0. Analysis using just the position of 
the first peak, also favors a ~ 34 . However, a complete analysis using the CMB data, in the sense of 
the analysis here developped for the H(z) data, presents new challenges. For example, in the references 
[321 133) , the parameter a has been restricted to positive values and only adiabatic perturbations were 
considered. Since the full perturbative analysis, at least at linear level, is necessary in order to obtain the 
CMB restrictions to GCG model parameters, some aspects already mentioned above point towards the 
necessity of a deeper analysis: the perturbative estimations for a show another peak of probability for 
very large a, either for matter power spectrum [20] and the Integrated Sachs Wolf effect [29], with this 
new peak located around a ~ 200; moreover, at perturbative level, the extension of the estimations to 
negative values of a is problematic, due to problems concerning imaginary sound speed. Some proposals 
to surmount this last difficult have been made in reference [13l E] but no clear scenario allowing to 
consider negative values of a has emerged until now. Then even if there is an overall agreement of 
the estimations made here and previous estimations obtained using CMB, the restrictions on the model 
parameters are not the same, and a new complete extensive analysis is required. Hence, the addition of 

1 The widening of the errors when the H(z) and SNIa data are crossed may indicate a tension between the data even 
taking into account the priors used. If there is this tension, it is small, and perhaps less important that other controversial 
aspects of the SNIa data, like the calibration issues, for example, those coming from the use of SALT or MLCS calibration 
methods 
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other cosmological probes does not provide a substantial change in the final parameter estimation. 




4 Conclusions 

The GCG model has four free parameters, when the flatness of the spatial section is imposed. Leaving 
the four parameters free, the unification scenario for dark matter and dark energy is not favored: the 
one-dimensional probability distribution for a extend to arbitrary negative values (if a m i n — —10), and at 
same time the prevision for n, m o becomes similar to the ACDM model. In fact, the crossing of these results 
with the ones comming from the perturbations (the matter power spectrum) [20) . allows to conclude: 
for any Chaplygin-based cosmology if the parameter f2 m o is free to vary the unification scenario, with 
&dmO ~ 0, is not favored. However, this problem can be alleviated if the allowed range of the parameter a 
is restricted to a > 0. With this choice, the unification scenario becomes the preferred one (see figure 3). 
Remark that this prior is obligatory for the matter power spectrum analysis in order to have real values 
for the sound speed. Again for the H (z) test, fixing Sl m o = 0.25, we also find a plateau for negative values 
of a. This result can be understood by remembering that large negative values of a implies a model very 
similar to the ACDM one. The imposition from the beginning of the unification prior (fi m o = 0.04) leads 
to a more consistent scenario, with a peak of probability for a near a = for any choice of ot m i n - Note 
that a joint analysis using also Supernovae as cosmological probes does not change the main predictions 
for this case. This is consistent with the corresponding power spectrum analysis, which implies also a 
peak of the probability around a = 0. In this sense, in order to allow the parameter a to assume negative 
values the prior of the unification scenario seems to be the only consistent one. 

The bulk viscous model appears as a second candidate for the unification scenario. Its background 
dynamics is similar to the GCG one. However, unlike the GCG, the perturbative analysis of the bulk 
viscous fluid supports the unification scenario. Having in mind the similarities between the bulk viscous 
model and the GCG one (section two), the constraints obtained in figure [5] agrees with the results for the 
GCG analysed in case (i). Due to the fact that the bulk viscous fluid has less pathologies or restrictions 
than the GCG one it seems put it as a more natural candidate for the unification scheme. One difficult 
with such model (based on the Eckart's formula ([2])) appears when dealing with the Integrated Sachs- 
Wolfe effect [17]. However, the viscous pressure has an alternative description in terms of the second 
order corrections to the equilibrium, namely the Mueller-Israel-Stewart theory [35] - This approach is 
very promissing for further investigation. 
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